The distance between Inherent Structures and the inluence of saddles on approaching 
the mode coupling transition in a simple glass former 
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We analyze through molecular dynamics simulations of a Lennard- Jones (LJ) binary mixture the 
statistics of the distances between inherent structures sampled at temperatures above the mode 
coupling transition temperature Tmct- After equilibrating at T > Tmct we take equilibrated 
configurations and randomly perturb the coordinates of a given number of particles. After that 
we find the nearby inherent structures (IS) of both the original and perturbed configurations and 
evaluate the distance between them. This distance presents an inflection point at Tu ~ 1 with a 
strong decrease below this temperature which goes to a small but nonzero value on approaching 
Tmct- In the low temperature region we study the statistics of events which give zero distance, 
i.e. dominated by minima, and find evidence that the number of saddles decreses exponentially 
near Tmct- This implies that saddles continue to exist even at T < Tmct- As at Tmct the 
diffusivity goes to zero, our results imply that there are saddles associated with nondiffusive events 
at T < Tmct- 
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What is the influence of the potential energy landscape 
on the dynamical properties of glass formers? Under- 
standing of this relation is crucial for the development of 
a comprehensible thermodynamics of supercooled liquids 
and glasses. In recent years much progress have been 
done in this direction (l). Simulating simple models of 
glass formers Sastry et al. |2) found evidence of four re- 
gions where the influence of the landscape is qualitatively 
different: (a) the high temperature region where the sys- 
tem is essentially a simple liquid, with exponentially de- 
creasing time correlation functions and free diffusion, (b) 
a landscape influenced region characterized by the onset 
of nonexponential relaxation for Tmct < T < Tu, (c) a 
region which Sastry et al. called landscape dominated in 
which activated events should be dominant and (d) the 
region for T < T g corresponding to the glass phase. In 
this work we will be concerned with the first two regions 
and the crossover between the second and the third one. 
In regions (a) and (d) the landscape seen by the system 
is essentially flat, of course for very different reasons: in 
(a) the energy of instantaneous configurations is much 
higher than typical potential energies and in (d) the sys- 
tem is able only to explore a strictly local region mainly 
through vibrations and diffusion is avoided. In (b) and 
(c) the ruggedness of the landscape has a strong influence 
on the dynamics and makes a thermodynamic description 
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very difficult to rationalyze. More recently a more quan- 
titative description of the role of the landscape has begun 
to emerge || g, 0, g|. An interesting way of look- 
ing at changes in the landscape topology in the different 
regions is to calculate the statistics of saddle points. By 
looking at the order of the saddles, i.e. the number of 
unstable directions n s , Angelani et al. Q found that this 
is a well defined function of temperature and found evi- 
dence that it goes to zero at a temperature which seems 
to coincide with the mode coupling transition tempera- 
ture Tmct [§• The immediate consequence of this result 
is that the mode coupling transition can be interpreted 
in terms of a topological change in the landscape: Tmct 
should mark the transition between a region dominated 
by saddles to a region dominated by minima. That this is 
qualitatively true is now reasonably well established but 
the quantitative identification of Tmct with the point 
where n s goes to zero needs to be independently verified 
for different systems and if possible by looking at several 
different quantities. One of these is the mean squared 
distance (MSD). In reference || the MSD was calculated 
between an instantaneous configuration and the nearby 
IS. The main conclusion was that at very low tempera- 
tures, tipically below Tmct, it goes linearly to zero as 
T — > 0, an evidence that the MSD comes essentially 
from harmonic vibrations. In reference Q it was shown 
the MSD between saddles and the corresponding minima 
suggesting a linear relation also between this MSD and 
the order of the saddles n s . 

Here we look at the distances between near inherent 
structures which give additional information on the re- 
lation between landscape topology and the dynamics of 
the system. We performed molecular dynamics Simula- 
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tions on a well known Lennard-Jones binary mixture of 
80 % particles of type A and 20 % particles of type B 
with e AA = 1.0, e BB = 0.5, e AB = 1.5, ctaa = 1.0, 
ubb = 0.88, ex ab = 0.8 at a density of 1.204 that has 
been extensively studied for 1000 particles by Kob et al. 
||Io| , [ll| who obtained Tmct — 0.435. In the present work 
we used N — N A + Nb — 250. We first equilibrate the 
system at temperatures T — 5,4,3,2,1,0.8,0.6,0.55,0.5 
and 0.47. We study 32 samples at each temperature fol- 
lowing different thermal treatments in order to have the 
same statistics than in reference [l(] and verify that ex- 
trapolating our data to lower temperatures the diffusivity 
should go to zero at the same value of Tmct- The dis- 
tance between IS is calculated as follows: once the system 
is equilibrated at a temperature T a typical configuration 
is taken. Then we make a damage on the equilibrated 
configuration by changing randomly the coordinates of rid 
particles chosen also randomly and the IS corresponding 
to the original configuration and the damaged one are ob- 
tained via conjugate gradients minimization. The MSD 
between both IS is defined as: 

D 2 (T,d) = (j2(rl S °-rl S «y) (1) 

i 

in which T is the temperature at which the system has 
been equilibrated, d is the amount of damage done to 
the equilibrated configuration, i.e. the number of parti- 
cles perturbed and the dislocation per particle and the 
braces mean an average taken over a number of indepen- 
dent equilibrium configurations and by making different 
damages on each configuration chosen, for example, by 
taking different particles to move for each configuration. 
As we want to explore the distance between near IS, 
i.e. those accesible dynamically at each temperature, the 
damage should be small, so we displace each one of the 
rid particles an amount |<5i*i | = 0.01 or \5vi\ = 0.10, both 
much smaller than the typical interparticle distance of 
the system p^] . Note that the distance D 2 goes to zero 
when the inherent structures of the original configuration 
and the damaged ones are the same, or in other words, 
when they are in the basin of a minimum, with no double 
well saddles between them. We found that at low tem- 
peratures an important fraction /o of the events gives 
zero distance, and so in order to get sensible results we 
normalize D 2 with the fraction of non zero events 1 — /o- 
Doing this, even in the case where only one event gives 
a non zero distance, the resulting D 2 will be finite and 
equal to this distance. Concerning the number of parti- 
cles to be damaged rid, we have checked that at low tem- 
peratures increasing rid from 2 to 40 only improves the 
statistics increasing 1 — /o but not changing qualitatively 

the behaviour of x _ J ' . So we take n<j=40 and keep it 
fixed in the following study, the amount of damage "d" 
in expression (1) can therefore be fully caracterized by 
the size of the dislocation \Sr\. For temperatures T < 0.6 
we have taken averages over 1000 configurations for each 
one of the 32 samples covering a time range greater than 
the relaxation time at each T. Then for each configura- 



tion we have also averaged over 4 and 8 different particles 
damaged for d = 0.10 and d = 0.01 respectively. 

A natural question is: is there some temperature at 
which D 2 (T,d) goes to zero? In this case this should 
coincide with the point at which the order of the saddles 
goes to zero as discussed above. 
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FIG. 1: Mean squared distance between inherent structures 
as function of temperature for two different damages d=0.1 
and 0.01. Inset: a zoom of the low temperature region. 

In figure [l] it is shown the normalized MSD as a func- 
tion of temperature for two different damages. As ex- 
pected the distance is very large at high temperatures 
and decreses upon decreasing T. Between T — 0.8 and 
T = 1 an inflection point is clearly seen below which the 
MSD decreases very quickly. This temperature range 
is the same at which Sastry et al. [§] found the onset 
of non-exponential behaviour in the relaxation functions 
signalling the begining of what they called the landscape 
influenced region. Our results are an independent confir- 
mation of the presence of that characteristic point but, 
more importantly, it gives a purely topological determi- 
nation of it. From the inset it is apparent that the dis- 
tance does not go to zero on approaching Tmct- 

In order to get a better insight into the landscape prop- 
erties of the low temperature region we show in figure || 
the probability of non zero events as a function of tem- 
perature in log-linear scale. It is evident the onset of a 
strong decrease in the non zero events at around T = 1. 
The behaviour of 1 — /o near Tmct is obviously depen- 
dent on the initial damage. But even for the smallest 
damage d = 0.01 (two orders of magnitude smaller than 
the typical interatomic distance) the number of non zero 
events does not seem to go to zero at Tmct- In the inset 
of figure H we show 1 — fo in a reduced temperature range. 
The points below T ~ 0.6 are perfectly aligned indicating 
an exponential decay of the number of non zero events 
towards T M ct of the form 1 — / = exp((T — T )/e) with 
Tq depending on the initial damage. We did several inde- 
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pendent checks in order to confirm this behaviour. The 
direct implication of this result is that the number of sad- 
dles does not vanish at Tmct in agreement with recent 
work ||| |j| and in apparent contradiction with previous 
one |Q|7 As we have independently verified in our sim- 
ulations that upon extrapolation to lower temperatures 
the diffusion constant should go to zero at Tmct, this, 
together with the previous result implies that there ex- 
ist saddles at T < Tmct associated with nondiffusional 
events. 
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FIG. 2: The fraction of events with non-zero final distance 
as a function of temperature for two different damages in log- 
linear scale. Inset: a restricted temperature range in order to 
show clearly the alignment of points for T < 0.6 

One can get a more direct evidence of the relation 
between the landscape topology and dynamics by look- 
ing directly at the particle displacements. The fact that 
particles cannot diffuse in a region where there is still 
minima-connecting saddles available may explain recent 
results obtained by Biichner et al. on the slowing down 
in the dynamics of supercooled liquids on approaching 
Tmct fl]| ■ In reference [l3|, the authors show that near 
Tmct the system visits a kind of valley containing sev- 
eral IS where the system is trapped for some time jump- 
ing between the different IS and where the mobility is ex- 
tremely small. Motion of particles in these valleys should 
be strongly collective and probably single displacements 
won't go beyond one interatomic distance. In order to 
test these hypotesis we have computed several displace- 
ment distribution functions. 

In figure || we show the distribution of squared dis- 
placements, i.e. each value of D 2 contributing to this 
histogram is a sum over the 250 particles which corre- 
sponds to each term contributing to the mean squared 
displacement of equation (Q) . We see that at high tem- 
peratures (T > 1) the maximum is located at D 2 greater 
than 40 which is the number of particles perturbed. This 
indicates that an important fraction of the system suf- 




FIG. 3: The distribution of squared displacements contribut- 
ing to the mean of equation (hj) for the different temperatures 
studied for damage d = 0.10. Note the presence of two max- 
ima for T=0.8 and T=l. Inset: the peak of the distributions 
for very small displacements. 

fer displacements at least of the order of the interatomic 
distance. An interesting fact emerges exactly at T = 1 
and is also present at T = 0.8: a second maximum at 
very low values appears. This is another indication that 
T = 1 is a characteristic temperature for this system be- 
low which the dynamics is landscape influenced as Sastry 
et al. named this region. In this crossover region three 
kinds of processes can be detected: (i) processes in which 
a fraction of particles can still move by distances of the 
order of an interatomic distance, (ii) processes where the 
particles can only move by tiny amounts and there is no 
particle moving an interatomic distance and (Hi) a third 
group of processes where particles do not move at all. 
The last group contributes with a delta peak at the origin 
which we have eliminated from our plots. The growing 
confinement of the particles as the temperature is lowered 
from T = 1 is more dramatic than what can be inferred 
from the growing peak at low values of the distance. In 
fact for T < 0.8 there is only one peak at D 2 < 0.5 (see 
the inset of figure ^) . If we note again that each contri- 
bution to D 2 is a sum over the 250 particles, this implies 
that on average no single particle moves a distance of 
the order of an interatomic distance. Of course, from the 
presence of the tails in the distributions for larger values 
of the distance one cannot discard a priori the existence 
of some particles that move by such amounts. Never- 
theless it seems more likely that this contributions come 
essentially from collective motions in which an important 
number of the particles move very little. A stronger sup- 
port for this scenario is given by the distribution of the 
largest displacements shown in figure ||. 

This figure gives also very interesting information. 
Note that for T > 1 the particle that moves more goes to 
a distance typically longer than one interatomic distance 
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PIG. 4: The distribution of largest displacements for several 
temperatures for damage d = 0.10. Dm ax is the displace- 
ment of the particle that experiences the largest displacement 
of all the particles in the system in a given process. Note the 
large dispersion in the values at T = 0.8 and T = 1. 

with very few contributions with Dm ax < 1- The vari- 
ance of the distributions grow and the height of the peak 
goes through a minimum as the temperature crosses the 
region where the landscape begins to influence the dy- 
namics. The large variance signals the appearance of an 
important fraction of particles which, although are the 
ones that move more than any other, have displacements 
confined to Dm ax < 1- Note that the distances between 
IS (figure [j]) present an inflection point between T = 0.8 
and T = 1 indicating a rapid decrease in the number of 
saddles available for diffusing. For T < 0.8 the maximum 
displacements are peaked around very small values less 
than one interatomic distance. 

Finally we asked the question: is there some order in 
the pattern of displacements per particle as a function 
of temperature? A partial but important answer can be 
done by ordering the N particles by decreasing value of 
their displacements. In figure |j| we show that in fact the 
patterns of displacements are similar for all temperatures. 

The first point to the left corresponds, for each tem- 
perature, to the mean over the particles with the largest 
displacements (first moment of distributions of figure [|) , 
the second point is the mean over particles with the sec- 
ond largest displacements, and so on up to the parti- 
cles with the smallest displacement. From this figure 
it is very clear that the overwhelming majority of par- 
ticles move in average less than 1 while only very few 
move more than one interatomic distance at high tem- 
peratures T > 1. For T < 0.8 all particles in average 
move very little. Interestingly, the average distance trav- 
eled by the particles with the largest displacement goes 
to zero very slowly (almost linearly) as T — > suggest- 
ing again a continuous variation of landscape properties 
on crossing Tmct- Also note that most particles move 
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FIG. 5: The pattern of displacements for the 250 particles 
and for several temperatures 



by some amount which means that the motion is highly 
collective. By looking at the relative difference between 
the first and second largest displacements (not shown) we 
obtained a practically constant and low value for T > 1 
followed by a pronounced grow of this difference for T < 1 
which is another indication of the drastic confinement of 
the system in this range of temperatures where only one 
particle moves much more than the others (although in 
a scale less than one interatomic distance). 

Recently Grigera et al. computed the distribution of 
particle displacements between a minimum and a saddle- 
connected neighbour in a soft spheres system and found 
evidence that, at low temperatures, the typical process 
is one where a large amount of particles travel very small 
distances, less than half an interatomic distance, and very 
few particles move by the order of one interatomic dis- 
tance. These were interpreted as activated processes in 
a region where there were almost no saddles available. 
This picture is qualitatively different from the one that 
emerges from the present work for the binary LJ system. 
Nevertheless, further work is needed to elucidate if dif- 
fusion in binary LJ system occurs as a cumulative effect 
of the small collective displacements or it is due to very 
rare activated events that, if present, should be confined 
to the tails of our displacement distributions. 

In summary we analyzed the information brought by 
the distances between near inherent structures relevant 
for the dynamics of a supercooled liquid near the mode 
coupling transition temperature. In this way the contri- 
butions coming from vibrations are automatically filtered 
out. It is possible to obtain rather precise information on 
the processes involved in the evolution in phase space and 
also in real space. The functional form of the MSD as a 
function of temperature and the evolution of the distri- 
bution of displacements and of maximum displacements 
define clearly a characteristic region near T = 1 below 
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which the number of saddles decays rapidly and the par- 
ticles become strongly confined. We found also that the 
number of saddles is exponentially small on approaching 
the mode coupling transition temperature but does not 
go to zero. In the low temperature region, in connecting 
two neighbouring IS, most of the times all particles move 
by amounts much smaller than the typical interparticle 
distance. The scenario that emerges for the low temper- 
ature dynamics is one in which the landscape is formed 
by a kind of metabasins [|| H [l^] in which there are still 
many inherent structures connected by low lying saddles 
which may exist even for T < Tmct (see also recent 
work by Heuer et al. ]l5||). Individual particles move 
by tiny amounts within a metabasin in a highly collec- 
tive way. From our results on the binary LJ system the 



mode coupling transition temperature Tmct continues 
to be characterized by a dynamical singularity and the 
possible connection of the transition itself with a sharp 
change in some topological property of the potential en- 
ergy landscape remains to be elucidated. 
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